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The stability of the spin-glass phase against a magnetic field is studied in the three and four 
dimensional Edwards-Anderson Ising spin glasses. Effective couplings J e a and effective fields H e g 
associated with length scale L are measured by a numerical domain-wall renormalization group 
method. The results obtained by scaling analysis of the data strongly indicate the existence of a 
crossover length beyond which the spin-glass order is destroyed by field H. The crossover length well 
obeys a power law of H which diverges as H — > but remains finite for any non-zero H , implying 
that the spin-glass phase is absent even in an infinitesimal field. These results are well consistent 
with the droplet theory for short-range spin glasses. 

PACS numbers: 75.10.Nr, 75.40.Mg, 05.10.Ln 



In spite of extensive studies for more than two decades, 
a basic problem on the field-temperature phase diagram 
of the short-range Ising Spin Glass (SG) is still controver- 
sial. The mean-field theory predicts the existence of the 
SG phase in a magnetic field up to a certain strength at a 
temperature below T c , the critical temperature in a zero 
field Q. On the other hand, the droplet theory @,[1,[3, 
a phcnomcnological theory for short-range SGs, predicts 
the absence of the SG phase even in an infinitesimal field. 

In experiments, this issue has been addressed by the 
study of Ising SG FexMni-^TiOs. Although the pres- 
ence of the SG phase in field was first reported @, the 
subsequent work concluded its absence by careful analy- 
ses of its relaxation time Q. The same conclusion was 
recently drawn by Jdnsson et al Q. From a theoret- 
ical point of view, numerical studies of the Edwards- 
Anderson (EA) short-range SG model [f| have yielded 
rather conflicting results: some data support the pres- 
ence of the SG phase in field @, [nj, and others support 
its absence ll,12,lj|. However, a recent numerical work 
on the correlation length has shown the absence of the SG 
phase in three dimensions at a low field H = 0.05 J fl3| . 
where J is the standard deviation of couplings {Jij}- 
Since this field is much below a critical field H c m 0.65 J 
suggested by a previous study [Io| . the result strongly 
indicates the absence of the SG phase. However, there 
still remains the possibility that the SG phase is stable at 
lower fields. Furthermore the physical mechanism which 
destroys the SG state remains to be clarified. 

In the present work, we study the EA SG in field in 
both three and four dimensions by the same numerical 
Domain- Wall Renormalization-Group (DWRG) method 
that we have recently used in studying the fragility of 
the SG equilibrium state against small changes in tem- 
perature or in bonds {Jij} Q- We also show some 
corresponding results of the Migdal-Kadanoff (MK) SG 



model, in which the absence of the SG phase in field has 
already been shown [TBI ]. For this model we can easily 
access huge sizes such as L w 10 10 . This compensates 
the Monte Carlo (MC) results on the EA model within a 
limited range of length scales. Quite interestingly, both 
the models exhibit the same scaling behavior. All the 
observables are scaled as a function of L/£ CT (H), where 
L is the system size and £ cr (H) a crossover length be- 
yond which the SG order is destroyed by field H. The 
existence of such a length scale is indeed predicted by the 
droplet theory. The observed scaling behavior indicates 
the absence of the SG phase even in an infinitesimal field. 

The droplet theory — Let us begin with a brief survey 
of the droplet theory 0, 0, 0] . Droplets are spin clusters 
which can be flipped with a low excitation energy. The 
typical excitation energy of droplets with length scale £ 
is assumed to scale as I , where 9 is the so-called stiff- 
ness exponent. Now let us consider to apply field H. 
Although we consider a uniform field for simplicity, the 
following argument is also valid for random fields. In 
Ising SGs, the spins in a droplet are either +1 or — 1 
with equal probability, implying that the order of Zee- 
man energy of droplets with size £ is Since the 
droplet theory provides some arguments which support 
the inequality 4=1 > 9 0, the theory claims that there 
exists a characteristic length, given as £ cr (H) ~ H^ 1 ^ 
with C, = d/2 — 9, that droplets larger than £ CI (H) are 
forced to flip by the field. As a result, the SG state at 
H = becomes unstable beyond £ cr (H). We call £ CI and 
£ the (field) crossover length and the crossover exponent, 
respectively. Furthermore the droplet theory claims that 
the system is paramagnetic beyond £ cr (H) as it happens 
in random field Ising model. Because £ C r{H) diverges as 
H — > but remains finite for any non-zero H, the droplet 
theory asserts the absence of the SG phase in any field. 

DWRG method and observables — Let us first explain 
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models for the DWRG method [13, El, E3] ■ Figure ED (a) 
shows the way to construct the hierarchical lattice for the 
MK SG. The lattice is made iteratively by replacing each 
bond with b d ~ 1 paths, where d is the dimension of the 
lattice. Each path consists of b bonds, and new (b — 1) 
spins (full circles) are inserted in between. We hereafter 
call the two outermost spins (Sl and Sr) boundary spins. 
The size L of the lattice is multiplied by b as the replace- 
ment is done once. Figure Q](b) is the lattice for the EA 
SG. The lattice is the same as that in ref. It consists 
of two boundary spins and L d spins on a rf-dimcnsional 
hyper-cubic lattice. The boundary condition is open in 
the direction along which the boundary spins lie, and is 
periodic in other directions. The Hamiltonian is 

H = -J2 I JijSiSj -^HiSi, (1) 

where the first term is exchange energies between two 
nearest neighboring spins and the second term is Zeeman 
energies by field Hi. 

In the DWRG method, we measure the effective cou- 
pling J G fj and the effective fields H c h defined by 

Z(S^Sn) = Tr'c- ff < s }/ T oce^( s - s *)/ T , (2) 



W e ff(S , L, Sr.) — — JcsSlSr — H^Sl — H^g' Sr, (3) 

where the trace in Eq. @ is for all the spins except Sl 
and Sr. In the MK SG, Z(Sl, Sr) is estimated exactly 
by taking the trace sequentially from the later generated 
spins to the earlier generated ones. For a detailed de- 
scription of the procedure, we refer the reader to [ID]. In 
the EA SG, on the other hand, probability P(5l,Sr) is 
measured by MC simulations [14| . Since J e e and the free- 
energy difference SF caused by twisting the two boundary 
spins are related by J c g = —SF/2 in zero field, we con- 
sider that J G ff represents the strength of the SG order. 
Because J e s is either positive or negative, we calculate 
the standard deviation of sample-to-sample fluctuations 
of effective couplings, <jj(L, H). We also measure that of 
effective fields a-n{L,H). The correlation between effec- 
tive couplings in zero field and those in field H is also 
estimated by the correlation coefficient 



r(R) , 



C(L,H) 



J Gff (L,g)J cff (L,0) 
a 3 (L,H)a 3 (L,0) : 



(4) 
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where 7TT is the sample average. Here J e ff(L,H) and 
J e ff(L, 0) are calculated for the same realization of {Jij}. 

Results in the MK SG — Figure [2] shows the size de- 
pendences of <7j and on in the four-dimensional MK SG. 
The values of couplings {Jij} are taken from a Gaussian 
distribution of mean and width 1. We apply random 
fields Hi of strength H by following the way in ref. [l5| . 
The pool method [l8| is used to access huge sizes such 
as L w 10 10 . The Zeeman energy (oh oc HL d / 2 ) over- 
whelms the effective coupling (<7j oc L e ) around L m 
H^ 1 ^, i.e., around the crossover length i CT (H) in the 
droplet theory. After the crossing, ctj exhibits roughly 
exponential decay and the exponent of cth changes from 
d/2 to (d — l)/2. These observations naturally lead us to 
the idea that H e ^aj, H s /^o~yi and C are scaled as func- 
tions of X = LH 1 ^ which we call the scaling variable. 
The idea is tested in Fig. [3J The data with different H 
and L nicely collapse into scaling curves. These results 
clearly show the existence of the crossover length £ cr (H). 

Results in the EA SG — Let us first explain some de- 
tails of our simulation. The values of Jij are either + J 
or — J with equal weights. Uniform field H is applied to 
all the spins except Sl and Sr. The number of samples 
is more than 1,500 for all the sets of (L,H). We use 
the exchange MC method to accelerate the equili- 
bration, and the method in ref. [l7| to overcome a hard- 
relaxing problem of the boundary spins which is origi- 
nated from their high connectivities. The temperature 
ranges used for the exchange MC are 0.5 J < T < 4.0J 
(T c w 1.1J Jp) for d = 3 and 1.0J < T < 4.5J 
(T c w 2.0J |22j |) for d = 4. We hereafter focus on the 
data at the lowest temperature which is well below T c . 
We set the MC step for thcrmalization and that for mea- 
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d=4, T=0.2 (T c ^1.6), H=10" 7 




FIG. 1: (a) The construction of the hierarchical lattice (b = 2, 
d = 3). (b) The lattice for the DWRG method in the EA SG. 



FIG. 2: Size dependences of <tj and an in the four dimensional 
MK SG at T w 0.125T C , where T c is the critical temperature 
in zero field. Strength of random fields is 1CP 7 . (A) Size 
dependence of oj for smaller L. (B) Size dependence of <th 
for smaller L. (C) The size where the crossover between oj 
and an occurs. (D) The value of <tj and <th at the crossing 
point. (E) Size dependence of <th for larger L. 
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surcmcnt to be equal. They are sufficiently (at least 5 
times) larger than the ergodic time to ensure the equili- 
bration, where the ergodic time is the average MC step 
for a specific replica to move from the lowest to the high- 
est temperature and return to the lowest one. As done in 
ref. 17| , we have also checked that MC runs starting from 
parallel and anti-parallel boundary spins yield identical 
results within error bars. 

Figures [4] and [5] show the results for d — 3 and those 
for d = 4, respectively. Since C is a dimensionless quan- 
tity, it is a function of only X = LH 1 ^ . We therefore 
estimate ( by fitting the data of C. By assuming the 
scaling relation ( = d/2 — 9 predicted by the droplet the- 
ory, 9 is estimated to be 0.29(3) for d — 3 and 0.60(2) 
for d = 4. Since t hey are close to both recent estima- 
tions by Boettcher [23| (0.24(1) for d = 3 and 0.61(2) for 
d = A) and direct estimations by linear least-square fits 
of ln[<7j(L, H = 0)] against ln(L) (0.28(3) for d = 3 and 
0.69(3) for d = 4), our data support the scaling relation. 

We next examine scaling properties of erj and cth by 
using the values of 9 determined above. The scaling rea- 
sonably works except crj for d = 4. The deviation sug- 
gests that the fields investigated are too large or/and the 
sizes are too small so that corrections to the scaling is not 
negligible. In fact, if we closely observe the scaling plot 
of er j for d = 3, we notice that the data with large H, say 
H > 0.48, systematically deviate from the master curve, 
while the data with small field arc scaled quite well. If 
we estimate 9 by forcing all the data to be scaled approx- 
imately, we get apparently larger values of 9 in both oj 



and oh- For example, 9 estimated from erj in such way 
is 0.36 ± 0.02 for d = 3 and 0.84 ± 0.04 for d = 4. 

The crossover behavior in the MK SG shown in Fig. [3] 
look very sharp in comparison with that in the EA SG. 
This is simply because the ranges of H and L, and so X , 
are quite wide. The enlargement of the crossover region 
such as to 0.1 < X < 10 is realized by choosing the pa- 
rameter ranges, for example, 0.015J < H < 0.95J and 
L = 3™ with 1 < n < 4. These ranges are close to those 
we are forced to choose for the EA SG. We then obtain 
quite similar results (not shown) as those obtained in the 
EA SG (Fig. |4]) . Within these ranges oj increases mono- 
tonically for smaller H and it decreases monotonically 
for larger H. Still, such data are reasonably scaled with 
the same value of 9 as indicated in Fig. [3] These observa- 
tions imply that the observed scaling behaviors in both 
the MK and EA SGs are essentially the same. 

Interpretation of results — We first consider the 
crossover behavior in oh- Since the field is not applied 
to the boundary spins, the effective fields and H^' 
originate only through the influence of the field applied 
to the bulk spins {Si}. For example, if the correlation 
between Si and Sl is positive, an upward field to Si tends 
to direct Si, upwards. Since the correlation can be either 
positive or negative depending on the site, the contribu- 
tion to the effective fields can also be either positive or 
negative. When the effective coupling J c g exists, Sl and 
Sr receive such random contributions, whose amplitude 
is proportional to H, from all the bulk spins. This yields 




1 10 d 1(f 10° 10° 1 




FIG. 3: Raw data (left panel) and their scaling plots (right 
panel) in the MK SG. The dimensions and the temperature 
are the same as those in Fig. [2] The value 8 = 0.77 obtained 
in ref. [l^] is used for the scaling. The two straight lines in 
the right-bottom panel are proportional to x d ^ 2 and x^ -1 ^ 2 , 
respectively. 



FIG. 4: Raw data (left panel) and their scaling plots (right 
panel) in the EA SG at T = 0.5J for d = 3. The value of 9 is 
estimated by fitting the data of C, and the value is used for 
the scaling plot of aj and that of oh. The slope of the line in 
the right-middle panel is obtained by linear least-square fits 
of ln[aj(L, H = 0)] against ln(L). 
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<7h which is proportional to HL d / 2 . When J e g vanishes, 
on the other hand, cth °c L^ -1 )/ 2 because the boundary 
spins, which interact with spins, receive contribu- 

tions only from spins around their surfaces. 

Now the meaning of our results are clear. As shown 
in Fig. [21 the Zeeman energy (cth) begins to overwhelm 
the effective coupling (oj) around L w £ cr (H) since an 
increases faster than aj. After that, erj rapidly drops to 
zero. This means that the SG order whose length scale 
is longer than £ cr (H) is destroyed by field. As described 
above, oh is kept increasing but with the change of its 
exponent from d/2 to (d — l)/2. The decay of C(L,H) 
to zero for larger scaling variable X indicates the vanish- 
ing of the correlation between the state in zero field and 
that in field. The observed scaling behavior consistently 
implies the absence of the SG phase in any non-zero field. 

Discussion and conclusions — Now let us argue on the 
fragility of the SG state against other perturbations. Ac- 
cording to the droplet theory, a change in the tempera- 
ture (or bonds) by the amount A gives rise to a new SG 
state which is decorrelated to the original one beyond 
the length scale £ c h(A). Here £ c h(A) is proportional to 
A" 1 ^* with (* = d s /2 - 9, d s being the fractal dimen- 
sion of droplets. This type of the fragility of the SG state 



is called temperature or bond chaos [24|. Our recent 
DWRG study [3] in the EA SG has indeed revealed the 
existence of ^ c h(A). A key observable in such studies is 
the correlation coefficient C of Eq. ((4]) with H replaced by 
A. As shown in Figs. 3, 4 and 5, the system exhibits sim- 
ilar fragility against the field perturbation. However, it 
must be noted that vanishing of erj for larger X strongly 
indicates that the magnetic field destroys the SG order 
itself in sharp contrast to the cases of temperature and 
bond perturbations. 




To conclude, we have studied the EA SG in field by a 
numerical DWRG method. The thermodynamic obscrv- 
ables of the EA SG are confirmed to follow the scaling in 
terms of the crossover length £ CX (H) as predicted by the 
droplet theory, the consequence of which is the absence of 
the SG phase in field. It should be noted that above the 
upper critical dimension d u = 6 different scenario may 
hold [25[ . We consider that all of our results concerning 
the temperature, bond and field perturbations provide us 
a strong evidence for the appropriateness of the droplet 
theory for the description of SGs in low dimensions. 
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FIG. 5: The same plot as that in Fig. U for the four dimen- 
sional EA SG at T = 1.0 J. 
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